Priedel oscillations induced by non— magnetic impurities 
in the two— dimensional Hubbard model 
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We study the interplay of correlations and disorder using an unrestricted Slave-Boson technique 
in real space. Within the saddle-point approximation, we find Friedel oscillations of the charge 
density in the vicinity of a nonmagnetic impurity, in agreement with numerical simulations. The 
corresponding amplitudes are suppressed by repulsive interactions, while attractive correlations lead 
to a charge-density-wave enhancement. In addition, we investigate the spatial dependence of the 
local magnetic moment and the formation of a magnetic state at the impurity site. 

PACS numbers: 71.27.+a, 71.55.-i 



I. INTRODUCTION 

The combined effect of correlations and disorder, espe- 
cially relevant in the context of high-Tc superconduc- 
tors, poses the theoretical problem of consistently in- 
cluding both mechanisms in a formal approacho The 
combined influence of disorder and many-body interac- 
tions can cause interesting new physical processes, and 
leads to unusual experimental properties. In addition 
to the drastic reduction of the transition temperature in 
cuprate superconductors, non-magnetic impurities can, 
for example, create local magnetic momentscTbJ vifhich 
extend over up to several lattice spacings and enrich the 
phase diagram by introducing new disordered magnetic 
regionsta. 

Electron density (Friedel) oscillations as a response to 
a local perturbation, i. e. the impurity, provides a clue 
to understanding the electronic properties (Fermi surface 
etc.) of a many-electron system. The period of the oscil- 
lation is determined by the Fermi surface, which may be 
defined by the singular points in the momentum distri- 
bution function. Especially in connection with the Fermi 
surface properties of high-temperature superconducting 
compoundal3 such as Bi2Sr2CaCu208+5, the Fermi sur- 
face evolution as a function of doping from the under- 
doped to maximally and over-doped regimes is an ex- 
tremely important, yet unresolved issueEJ. A systematic 
study of Friedel oscillations, both in experiment and in 
the generic 2-D Hubbard model for the high-Tc cuprates, 
could, in principle, provide the required information on 
the Fermi surface and how it changes from the extremely- 
low doped case (with the possibility of "hole pockets", i. 
e. "small" Fermi surfaces) to the larger doped situation 
with a large Fermi surface, with the latter being in ac- 
cordance with Luttinger's-theorem. 

In a previous publicationlij, a diagrammatic scheme was 
proposed that uses the T-matrix formulation to derive 
an effective scattering potential which includes the ef- 
fect of correlations as well as the bare scattering poten- 
tial. As the first step of a systematic approach, we focus 
here on the one-particle effects of potential scattering in 



correlated systems by formulating a static ansatz using 
an unrestricted Slave-Boson (USB) method, which cor- 
responds to the effective-medium renormalization of the 
scattering potential. 

The calculations presented here are done for the slightly 
overdoped case, i.e. 17% doping away from half fiUing (i. 
e. the insulating case for the 2-D Hubbard model). In 
this case, we observe the 2fci?-oscillations, essentially un- 
renormalized from the ones of the non-interacting case. 
Not only periodicity, but also the ampHtude of the Friedel 
oscillations are shown here to be in excellent agreement 
with Quantum Monte Carlo (QMC) simulations. Cor- 
responding calculations and comparisons with QMC for 
smaller doping are presently carried out, in order to test 
the simple effective medium approach also in the most 
interesting doping regime. 



II. FORMALISM 

We start with the Slave-Boson technique as formulated 
by Kotliar and RuckcnsteinEJ. An elementary introduc- 
tion to the Slave-Boson formalism for disorder-free appli- 
cations can be fojind, for example, in the review articles 
by Arrigoni et alnA and Dieterich et alcn. The canonical 
transformation 



(1) 



of the fermionic annihilation (creation) operator cijj in- 
troduces new bosonic operators e\^\ p[^J , d'l'^ on a site 
i, which can be empty, singly occupied, or doubly occu- 
pied. Within the physical subspace the Hamiltonian can 
be written as 



SB = E 4/1/,.%. + UJ2 4d^ + VoJ2 /o./o.- 



(2) 

Here the impurity is modeled by a one-particle potential 
Vb acting at the site r = of the lattice. This model can 
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be used, for example, to-describe non-magnetic scatter- 
ers in high-Tc material^. 

In the paramagnetic (spin-independent) saddle-point ap- 
proximation of the path-integral, one calculates the par- 
tition function Zstat = e""^^*"' by determing the station- 
ary value of the effective action 



Sstat = -Trlog [5,j{dr - M + Vo5,o + Af ^) -tz,z, 



^USB 



(8) 



+2p,(Af'-Al^0p. + e.Al^^eO, 
which is given by the saddle-point equations, i. e. 

^ = 0, 5.e{e.,p.,d.,A«,A,r}. 



i(2) 



(3) 



(4) 



The derivatives must be calculated with respect to all 
site-dependent bosonic variables and all Lagrange mul- 
tipliers Aj-^^^\ which are introduced in order to restrict 
the enlarged, mixed bosonic-fermionic Hilbert space r-to 
the physical subspace within the functional integralE3. 
Note that we use the standard renormalization factors 
= /, L — /, \ — 2 Eq. (P) to guarantee the 

correct U limitnJ. On a 10 x 10 lattice and working 
at fixed average particle density n, one obtains 85 self- 
consistency equations by exploiting the spatial symmetry 
of the bosonic variables which are directly connected to 
physical observables. Within the Slave-Boson formalism, 
the particle density and the local magnetic moment are 
given by 



and 



n{ri) = 2dj + 2pl 



m{fi) ^ 2p1 



(5) 



(6) 



Compared to mean-field Slave-Boson approaches, which 
can be restricted to homogeneous, spiral or bipartite con- 
figurations, it takes considerably more numerical effort 
to find the extreme value of Sstat in the high dimen- 
sional phase space encountered here. The strength of 
the Slave-Boson technique is that many-particle correla- 
tions are included via the z-factors in the effective masses 
of the quasi-free fermions. For the broken translational 
symmetcY considered here, this renormalized mean-field 
methodEj also takes into account site-dependent hopping 
amplitudes i^J^ = tZiZj. 

By examining the fermionic propagator G/ occurring in 
Eq. (0), one obtains the Dyson equation 



= - [{dr ~ M)<5y - t z*z,] - [VAo + Xf^)5,, 
- (j-p^ - 2. , 

where the self energy is defined by 



(7) 



In an effective-medium or Ha|i±i;ee-like description of the 
scattering process one findsEjO 



^Hartree 



(9) 



with site-dependent particle densities (ni). Using the 
non-interacting propagator G^^ — —{dr — n)Sij — t z*Zj, 
Eq. can be represented diagrammatically as in Fig. 
1 (compare also to Fig. 1 of Ref. 17). Therefore, the 
correspondence of the Hartree- and the USB-method can 
be seen in the additional one-particle potential due to 
the interplay of correlations and disorder on the static 
"mean-field" level considered here. 



III. RESULTS 

In the past, the Slave-Boson method has been used as 
a powerful tool to calculate local quantities, which have 
turned out to be. in excellent agreement with numeri- 
cally exact resultsOouJ. In order to check the reliabil- 
ity of the USB-method in the context of the extremely 
large system of coupled saddle-point equations, we com- 
pared our results to Quantum-Monte-Carlo calculations. 
The QMC-simulations suffer here from the problem of 
increasing numerical instability due to the broken trans- 
lational symmetry, which limits the system size and the 
lowest accessible temperatures. Fig. 2 shows the com- 
parison of the USB-data with QMC results for an 8 x 8- 
lattice. The charge density for U — At (in the following, 
we set t = 1) is plotted as a function of the distance r/a 
from the impurity site, where a is the lattice constant. 
The saddle-point solution is almost identical to the QMC 
result despite the fact that the latter calculation is carried 
out at a relatively high temperature {(3 — 6,T — 1/6). 
Fig. 3 displays the USB oscillations of the charge density 
in the vicinity of the impurity at an inverse temperature 
/3 = 100 for various values of the correlation strength 
U. The solid line represents the non-interacting sys- 
tem and exhibits the well-known 2A:ir-oscillations, cor- 
responding here to a waveleagth of A « 1.5a. While the 
2A:ir-periodicity is unaffectedEa, repulsive interactions di- 
minish the amplitude of the Friedel oscillations. They 
lead for strong correlations {U > 10) eventually to a 
short-range, monotonically decaying behavior. The rea- 
son for this is that for increasing U the renormalized 
(pseudo-fermionic) particles with correlation-enhanced 
effective masses become more and more ineffective in con- 
tributing to screening processes. For attractive interac- 
tions, the impurity locally breaks particle-hole symme- 
try and causes large Friedel amplitudes. This is due to 
the tendency of the attractive Hubbard model to form 
a charge-density-wave instability, iyhich occurs in the 
unperturbed system at half-fiUingEl'Ell. The oscillations 
present at r/a > 5 are due to interference effects at the 
corners of the 10 x 10 lattice. They are irrelevant for, for 
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example, averaged random impurity sites. 
Fig. 4 exhibits the influence of the potential strength Vq 
on the corresponding charge density behavior. Increas- 
ingly attractive impurities enhance the amplitude until 
the maximum density (n —^ 2) is reached at the im- 
purity site r = (inset of Fig. 4). Stronger potentials 
decrease the energy of the bound states below the quasi- 
particle band. A symmetric repulsive-attractive effect of 
the impurity can be found in the oscillations for Vq = ±1 
around the average density (n = 0.81). Further USB- 
calculations showed that higher temperatures cause a de- 
crease of Friedel oscillations, which is in accordance with 
the smearing out of the jump in the Fermi function at 
kp, which determines, fpt example, the integral bond in 
a continuum calculationE3 of n{r). This is also the origin 
of the moderate oscillations in Fig. 2 at /3 = 6. We note 
that no significant deviations from the above described 
behavior were found in the investigated metallic phase at 
0.7 < n < 0.9. 

We now concentrate on the one-particle spectra obtained 
from the USB-calculations. Using the pseudo-fermionic 
propagator, we obtain the density of states 

D{uj) = -^^liRGf{m,uj + ir]), 

This function, with the (5-functions replaced by rj- 
broadened functions, is plotted in Fig. 5. The broken 
translational symmetry lifts the degeneracy of various 
quasi-particle states of the unperturbed system (sohd 
line) and produces a single bound state (arrow in Fig. 
5) separated from the band. Note that this band has a 
smaller width than in the non-interacting case due to the 
Slave-Boson renormalization effects. The bound state 
possesses s-wave symmetry (see, for example, the first 
state in Fig. 6), as does the bare scattering potential 
Therefore, within our static potential renormalizationEj, 
the effective impurity potentiaO is not strong enough to 
cause the pjmviously proposed (and verified in numerical 
simulationsEj) anisotropic scattering processes, resulting, 
for example in bound states of higher symmetries. These 
extended bound states can be found by introducing an 
additional bare one-particle potential Vi at the neighbor- 
ing sites, which produces states of all the possible irre- 
ducible representationsEj. Fig. 6 shows the correspond- 
ing amplitudes of the pseudo-fermionic (bound state) 
wave-functions \^m) — J2i'^rnifj\vac), belonging to the 
six lowest lying energies Em- 

Another interesting feature, found in the USB-results, 
occurs in the magnetic structure influenced by the non- 
magnetic impurity. We first investigate the local mag- 
netic moment to(0) at the impurity site r = 0. In Fig. 
7, this quantity is plotted (for fixed U and average parti- 
cle density n) as a function of the potential strength Vq 
and the corresponding local particle density ?i(0). The 
maximum of m(0) at Vq ~ —1.0 corresponds to a locally 



half-filled impurity site, i. e. n(0) = 1, and signals the 
formation of a magnetic state. This situation is analo- 
gous to the mean-field description of the Single-Impurity 
Anderson model and was also found in exact diagonal- 
ization studies of a Hubbard chaiji. containing a single 
impurity by Hallberg and Balseirocj. In order to evalu- 
ate the size of the magnetic moment, one should perform 
spin-dependent USB-calculations, which could then also 
be extended.|tfl the symmetry-broken phases found near 
half-fiUingOLj. The paramagnetic subspace to which we 
restrict the solution appears sufficient for the questions 
investigated here. Substituting U + Uq for the local in- 
teraction strength f7 at r = 0, one shifts the magnetic 
transition point by a small Uq, which can be seen in the 
inset of Fig. 7. This means that the appearance of a 
magnetic state depends on a fine tuning of the system 
parameters, which could be one reason for the experi- 
mental disagreement coucerning the size and existence 
of a magnetic momentErHl. Strong potentials eventually 
lead to the formation of a = spin-singlet (in the 
limit of no or double occupancy), which corresponds in 
Fig. 7 to a minimal local moment at the impurity site. 
Interestingly, the oscillations of the charge density are 
also accompanied by spatial variations of the magnetic 
moment. Using the parameters relevant for Fig. 7, Fig. 8 
shows the magnetic structure in the vicinity of the impu- 
rity. For Vq = —1, the magnetic state at r = produces 
a relatively small amplitudes m(r). The stronger attrac- 
tive potential, Vq = —4, creates a singlet bound state 
at r = and leads to a maximal moment at a distance 
of r = V2. This can be connected to the experimen- 
tally detected spatially extended (para-)magnetic mo- 
ments around Zn-impurities. In a theoretical picture, 
this feature of the metallic phase can be seen as a rem- 
iniscence of correlation-induced bound states- carrying 
localized magnetic moments, at half-filling23^o. 

IV. CONCLUSION 

In summary, we have extended the Slave-Boson mean- 
field method to systems with a spatially inhomogeneous 
phase space. This technique can be used to calculate 
single-particle properties of correlated-lattice models 
which include a non-magnetic impurity. Formally, this 
method corresponds to a static renormalization of the im- 
purity potential in analogy to a diagrammatic effective- 
medium Hartree approach to impurity scatteringlj. 
We have solved numerically the complicated system of 
85 coupled saddle-point equations for a 10 x 10 Hubbard 
model in the metallic phase. The USB-results show, in 
accordance with QMC-data, Friedel oscillations of the 
charge density in the vicinity of the impurity. The am- 
plitude is reduced by repulsive correlations, while the at- 
tractive system display a charge-density-wave enhance- 
ment. The system parameters can be chosen, so that 
a magnetic state appears at the impurity site, signaled 
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by the maximal local magnetic moment. In addition, 
the range of the magnetic correlations around the impu- 
rity increases with increasing potential strength, being in 
qualitative accordance with experimental results. 
The above formalism and its results can be viewed as 
the first step of a systematic investigation of correlation 
effects on impurity scattering. Higher diagrammatic or- 
ders of maruii-body interactions will be studied in a future 
publicationcJ. 
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Figure captions 

FIG. 1: Diagrammatic representation of the self- 
consistent pseudo-fermionic one-particle propagator Gf 
(double line). The scattering by the impurity is depicted 
by the cross, the zig-zag line represents the Slave-Boson 
self-energy, i. e. the Lagrange parameter A*^^). 

FIG. 2: Radial modification of the charge density caused 
by a (5-potential Vq = —8.0 The data points of the 
Slave-Boson calculation fit the QMC-distribution for the 
Hubbard-model at f/ = 4, /3 = 6 and an average density 
n = 0.83. 

FIG. 3: Correlations effects on the low-temperature 
(/3 = 100) Priedel oscillations induced by a 5-potential 
of strength Vq = -4.0. 

FIG. 4: USB-calculations for different impurity poten- 
tials in the Hubbard model {U = A, (3 = 100). Repulsive 
(attractive) interactions reduce (enhance) the amplitude 
of Friedel oscillations. 

FIG. 5: USB-spectra of the pseudo-fermionic propaga- 
tor. The density of states shows an impurity-induced s- 
wave bound-state below the renormalized quasi-particle 
band. The broken translational symmetry creates states 
between the unperturbed (Vq = 0) single-particle energy 
levels. 

FIG. 6: Wave functions of the six lowest pseudo- 
fermionic eigenstates caused by an extended impurity 
(Vb = -8.0, V\ = -4.8). The wave functions exhibit 
Px/y~ (degenerated), s-, dx2-y2- and s-wave symme- 
try, in order from upper left to lower right. 

FIG. 7: Charge density and magnetic moment at the im- 
purity site (r = 0). The scattering potential acts in the 
(7 = 4 Hubbard model at (3 = 100. The maximum of 
m(0) and the half-filled site displays the formation of a 
magnetic state. 

FIG. 8: The local magnetic moment in the vicinity of 
the impurity. The spatial magnetic distribution and its 
maximum depends on the strength of the scattering po- 
tential Vq. 
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